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According to the no-hair theorem, astrophysical black holes are uniquely characterized by their 
masses and spins and are described by the Kerr metric. Several parametric spacetimes which deviate 
from the Kerr metric have been proposed in order to test this theorem with observations of black 
holes in both the electromagnetic and gravitational-wave spectra. Such metrics often contain naked 
singularities or closed timelike curves in the vicinity of the compact objects that can limit the 
applicability of the metrics to compact objects that do not spin rapidly, and generally admit only 
two constants of motion. The existence of a third constant, however, can facilitate the calculation 
of observables, because the equations of motion can be written in first-order form. In this paper, 

I design a Kerr-like black hole metric which is regular everywhere outside of the event horizon, 
possesses three independent constants of motion, and depends nonlinearly on four free functions that 
parameterize potential deviations from the Kerr metric. This metric is generally not a solution to the 
held equations of any particular gravity theory, but can be mapped to known four-dimensional black 
hole solutions of modihed theories of gravity for suitable choices of the deviation functions. I derive 
expressions for the energy, angular momentum, and epicyclic frequencies of a particle on a circular 
equatorial orbit around the black hole and compute the location of the innermost stable circular 
orbit. In addition, I write the metric in a Kerr-Schild-like form, which allows for a straightforward 
implementation of fully relativistic magnetohydrodynamic simulations of accretion hows in this 
metric. The properties of this metric make it a well-suited spacetime for strong-held tests of the 
no-hair theorem in the electromagnetic spectrum with black holes of arbitrary spin. 

PACS numbers: 04.50.Kd,04.70.-s 


I. INTRODUCTION 

According to the no-hair theorem, isolated and station¬ 
ary black holes in general relativity are uniquely charac¬ 
terized by their masses M and spins J and are described 
by the Kerr metric. This metric is the unique station¬ 
ary, axisymmetric, asymptotically flat, vacuum solution 
of the Einstein field equations which contains an event 
horizon but no closed timelike curves in the exterior do¬ 
main 0 , 1 - Thanks to the no-hair theorem, all astro- 
physical black holes are expected to be Kerr black holes. 
While observational evidence suggests the existence of 
event horizons in astrophysical black holes (see the dis¬ 
cussion in, e.g.. Ref. Q), a proof of the validity of the 
no-hair theorem is still lacking. 

Astrophysical black holes, however, will not be per¬ 
fectly stationary or exist in perfect vacuum because of the 
presence of other objects or helds such as stars, accretion 
disks, or dark matter, which could alter the Kerr nature 
of the black hole. Nonetheless, under the assumption 
that such perturbations are so small to be practically un¬ 
observable, one can argue that astrophysical black holes 
are indeed described by the Kerr metric. This is the as¬ 
sumption I make in this paper. This is usually a good 
approximation for supermassive black holes in the centers 
of galaxies and for stellar-mass black holes in X-ray bina¬ 
ries, which are typically separated from their respective 
companion stars by ~ 1 AU (see, e.g.. Ref. i)- 

Several model-independent strong-held tests of the no¬ 


hair theorem have been suggested using gravitational- 
wave observations of extreme mass-ratio inspirals (EM- 
RIs) fsl-lla a nd electromagnetic observations of accretion 
hows pill - 1^ . These tests are designed in a phenomeno¬ 
logical approach that encompasses large classes of mod¬ 
ihed theories of gravity instead of focusing on any par¬ 
ticular gravity theory. In this approach, the underlying 
theory is usually unknown, but it is assumed that par¬ 
ticles in this theory move along geodesics. The goal is, 
then, to gain insight in this theory through observations. 
See Ref. for reviews of such tests. 


Other tests of the no-hair theorem exist that include 
the observation of gravitational ringdown radiation of 
perturbed black holes after a merger with another object, 
which tests whether the end-state of the merger is a Kerr 
black hole [s^, as well as several weak-field tests in the 
electromagnetic spectrum such as those obtained from 
the monitoring of close stellar orbits around Sgr A* [sij 
and pulsar/black-hole binaries 32|. 

While it is sufficient for tests in the weak-field regime to 
rely on a parameterized post-Newtonian approach (PPN; 
e.g., (H), the mo del-independent strong-held tests de¬ 
scribed above require a modihed spacetime which devi¬ 
ates from the Kerr metric in parametric form. Several 
such parametric frameworks have been created, within 
which possible observational signatures of a Kerr-like 
black hole can be explored (e.g., i i 0,0,0IH). 
These metrics can deviate slightly to severely from the 
Kerr metric, and observables can be studied in terms of 
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one or more free parameters. All of these metrics reduce 
to the Kerr metric if the deviations vanish. Since gen¬ 
eral relativity has to date only been marginally tested 
in the strong-held regime [s^, deviations from the Kerr 
metric could be either small or large as long as they are 
consistent with current weak-held tests (see [ 13 ). 

If a deviation from the Kerr metric is detected, there 
are two possible interpretations. Within general relativ¬ 
ity, the object in this case cannot be a black hole, but 
is instead a stable stellar conhguration or else an exotic 
object [ 1 , However, if the compact object is known 
to possess an event horizon, then the no-hair theorem 
is falsihed and four-dimensional general relativity is only 
approximately valid in the strong-held regime. 

Because of the no-hair theorem, all parametrically de¬ 
formed Kerr spacetimes have to violate at least one of 
the assumptions of this theorem. These metrics are gen¬ 
erally stationary, axisymmetric, and asymptotically hat. 
If such a metric is also Ricci-hat it either does not har¬ 
bor a black hole or contains singularities or regions with 
closed timelike curves outside of the event horizon. Such 
pathologies may also exist, if the metric is not Ricci- 
hat. The presence of pathologies in a given metric ham¬ 
pers using the metric for strong-held tests of the no-hair 
theorem in both the gravitational-wave and electromag¬ 
netic spectra, because physical processes in the immedi¬ 
ate vicinity of the black hole at radii comparable to the 
innermost-stable circular orbit (ISCO) cannot be mod¬ 
eled properly. Future gravitational-wave detectors, how¬ 
ever, will be most sensitive to EMRIs that occur at radii 
roughly in the range between the innermost stable orbit 
and 10 — 20M and accretion disks are often modeled to 
terminate at the ISCO. 

For certain applications such as ray-tracing simulations 
of the electromagnetic radiation emitted from geometri¬ 
cally thin accretion hows a cutoff radius can be imposed 
which acts as an artihcial horizon and shields the ad¬ 
verse effects of the pathological regions from distant ob¬ 
servers. While, in principle, information on an enclosed 
pathology will be encoded in the physical boundary con¬ 
ditions of the cutoff which could still affect the causal 
past of particles or observers outside of the cutoff, in 
practice, any photon that reaches the cutoff radius can 
be excluded from the simulation without altering the ob¬ 
served signal (see the discussion in, e.g.. Ref. [2^). For 
rapidly-spinning black holes, however, the ISCO can be 
arbitrarily close to the event horizon. Therefore, such a 
metric can only be used for strong-held tests with at most 
moderately-spinning black holes as long as the ISCO lies 
outside of the cutoff radius. See Ref. for further dis¬ 
cussion. 

Nonetheless, it is possible to design Kerr-like metrics 
that are free of such pathologies outside of the central 
object [l3, [l^. These two metrics parameterize devia¬ 
tions from the Kerr metric in generic form. While the 
metric proposed in Ref. [lj[ generally describes a black 
hole, the metric of Ref. harbors a black hole only 
for small perturbations away from the Kerr metric and 


generally describes a naked singularity. In this case, how¬ 
ever, it is always possible to choose a cutoff radius be¬ 
tween the surface of the naked singularity and the ISCO, 
which still allows for the study of the electromagnetic 
emission from, e.g., geometrically thin accretion disks in 
this metric [s^. The metric of Ref. [3 also admits an 
approximately (i.e., for small deviations from the Kerr 
metric) conserved third constant of motion in addition 
to the two associated with stationarity and axisymme- 
try, which is useful for the construction of approximate 
EMRI waveforms [l^. These metrics, therefore, can be 
used for strong-held tests of the no-hair theorem even for 
large values of the spin [s^. 

The structure of these metrics, however, has several 
disadvantages. The metric of Ref. depends only on 
one set of deviation parameters. In general relativity, sta¬ 
tionary, axisymmetric, asymptotically hat, vacuum met¬ 
rics can be written in terms of four independent functions 
[40l| , and one would expect that such metrics in more gen¬ 
eral theories of gravity should depend on at least four in¬ 
dependent functions and, thus, on four sets of deviation 
parameters. Likewise, it not obvious how to implement 
three-dimensional fully relativistic hydrodynamic simu¬ 
lations of accretion hows in the presence of a cutoff [s^ . 
The metric of Ref. [l^ depends only linearly on four sets 
of deviation parameters. This may not be the optimal 
form to parameterize such deviations, because there is 
no reason to believe that they are small. 

Such a dependence on the deviation parameters is suf- 
hcient for tests of the no-hair theorem with EMRI obser¬ 
vations as long as the deviations are assumed to be small, 
because an extension to include large deviations would re¬ 
quire knowledge of the strong-held radiative dynamics in 
the underlying theory (e.g., wave emission and polariza¬ 
tion, radiation reaction, matter coupling). Strong-held 
tests of the no-hair theorem in the electromagnetic spec¬ 
trum, however, are not a priori limited to the study of 
small deviations from the Kerr metric, because these are 
performed in a stationary black hole spacetime, where 
the metric serves as a hxed background. Therefore, the 
dynamical properties of the gravity theory are not im¬ 
portant and the underlying held equations do not have 
to be known. 


Electromagnetic tests of the no-hair theorem, then, 
only require a suitable accretion how model, for which 
additional assumptions regarding the coupling of matter 
to electromagnetic helds in the modihed theory have to 
be made. The simplest approach is to assume that all 
electromagnetic interactions are governed by the laws of 
Maxwell electrodynamics. This is a reasonable choice, 
given that, at present, there is no evidence for violations 
of this theory at the classical level (see, e.g., [13 )■ In 
many geometrically thin accretion disk models, the disk 
plasma is constrained to move along circular equatorial 
orbits and photons are emitted by plasma particles in 
their respective local rest frames according to the pre¬ 
scriptions of elementary atomic physics without any sub¬ 
sequent interactions. The resulting observed spectra can 
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be modeled using a variety of different ray-tracing algo¬ 
rithms (e.g., El). 

In fully relativistic magnetohydrodynamic simulations 
of accretion flows, the dynamics of the flow plasma and 
electromagnetic fields are governed by the equations ex¬ 
pressing the conservation of particle number and the 
stress-energy of the particles. Such simulations have, 
so f ar, only been performed in general relativity (e.g., 
[4^ SjI), where the flow particles are usually described 
as a perfect fluid, but these can likewise be carried out 
in a suitable Kerr-like metric. In many numerical cal¬ 
culations involving either geometrically thin or thick ac¬ 
cretion flows, however, a linear dependence on the devi¬ 
ation parameters cannot be enforced, because the equa¬ 
tions that govern the evolution of particles and fields such 
as the geodesic equations are usually nonlinear. It is, 
therefore, desirable to generalize the metrics proposed in 

Refs. [13, il. 

The existence of three constants of motion makes a 
metric special, because it allows for the separation of 
the geodesic equations (the fourth integral of motion is 
trivially provided by the normalization condition of the 
4-momentum of the test particle). Otherwise, geodesic 
orbits are generally chaotic (see, e.g., Refs. [43,[45|). The 
Kerr metric possesses such a third constant of motion, the 
Carter constant, which is associated with a second-rank 
Killing tensor [^. Brink [l3, 113, 113 investigated the 
existence of a third constant of motion in general sta¬ 
tionary, axisymmetric, vacuum spacetimes and showed 
that for a certain subclass of these metrics a fourth-rank 
Killing tensor exists. In Newtonian gravity and Maxwell 
electrostatics, a Carter-like constant exists if the source is 
axisymmetric and only has even multipole moments that 
have the same structure as the multipole moments of the 
Kerr metric (c.f., Ref. [l^l and references in Ref. [43|1. 
The general-relativistic analogue of such a matter con¬ 
figuration, however, does not possess a Carter-like con¬ 
stant [43| . 

In this paper, I design a new Kerr-like black hole met¬ 
ric which suffers from no pathologies in the exterior do¬ 
main, admits three independent, exact constants of mo¬ 
tion, and depends nonlinearly on four independent de¬ 
viation functions that measure potential deviations from 
the Kerr metric in the strong-field regime. This metric 
is not a solution to the Einstein field equations or any 
concrete modified theory of gravity, but it can serve as a 
phenomenological framework for strong-field tests of the 
no-hair theorem in general classes of gravity theories. 

In order to design a metric with these properties, I 
carefully introduce arbitrary deviations from the con- 
travariant Kerr metric in a manner that preserves the 
separability of the Hamilton-Jacobi equations. From this 
ansatz, I obtain a covariant metric which describes a 
black hole and which is regular everywhere outside of 
the event horizon. I further simplify this metric by re¬ 
quiring that it is asymptotically flat, possesses the cor¬ 
rect Newtonian limit in the non-relativistic regime, and 
is consistent with all current weak-field tests in the PPN 


formalism. 

First I design the metric in Boyer-Lindquist-like coor¬ 
dinates and derive a Carter-like constant as well as so¬ 
lutions of the Hamilton-Jacobi equations. By choosing 
expansions of the deviation functions in suitable power 
series I write the metric in an explicit form and impose 
consistency with the current PPN constraints. I then 
show that the event horizon of this metric is identical to 
the Kerr event horizon. 

I proceed to analyze circular equatorial geodesic orbits 
in this metric. I derive expressions for the energy, an¬ 
gular momentum, and orbital frequencies of particles on 
such orbits as a function of the mass, spin, and devia¬ 
tion parameters of the black hole and show that they are 
significantly modified with respect to their forms in the 
Kerr metric. From the energy I calculate the location of 
the ISCO and show that its location is likewise shifted 
significantly. I further derive expressions of the principal 
null congruences in this metric and find a transformation 
of the Boyer-Lindquist-like coordinates to Kerr-Schild- 
like coordinates. With this transformation I remove the 
coordinate singularity at the event horizon which allows 
for a straightforward implementation of fully relativistic 
magnetohydrodynamic simulations of accretion flows in 
this metric. 

Finally, I relate this metric to known four-dimensional, 
analytic black hole solutions of modified theories of grav¬ 
ity. I map the new metric in a form appropriately lin¬ 
earized in the deviation parameters or the spin to the 
black hole solutions in Einstein-Dilaton-Gauss-Bonnet 
[soil and Chern-Simons gravity and as well as to the 
Kerr-like metric of Ref. [I4| written in one of its forms in 
Ref. [l3|. In the latter case, I find exact recursion rela¬ 
tions between the deviation parameters of both metrics, 
which I verify explicitly at the first few nonvanishing or¬ 
ders. I also point out that the black hole solution [s^] 
in Randall-Sundrum-type braneworld gravity [^ can be 
trivially included by choosing the Kerr-Newman metric 
as the starting point of my analysis. 

The paper is organized as follows: In Sec. m I design 
the metric. In Sec. mi I calculate the locations of the 
event horizon. Killing horizon, and ergosphere and deter¬ 
mine the conditions on the deviations functions so that 
the exterior domain is regular. I analyze geodesic motion 
in the equatorial plane and the location of the ISCO in 
Sec.HV] and derive the principal null directions and the 
Kerr-Schild-like form of the metric in Sec. El In Sec. I VII 
I map the metric to known metrics from the literature. 
I formulate my conclusions and discuss astrophysical ap¬ 
plications in Sec. lVHI Throughout, I use geometric units, 
where G = c = 1, unless I state it explicitly. 


II. DESIGN OF A BLACK HOLE METRIC 
WITH THREE CONSTANTS OF MOTION 

In this section, I design a new class of stationary, 
axisymmetric, asymptotically flat metrics that describe 
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spinning black holes and that admit three independent 
constants of motion. The metric elements depend on 
the mass and spin of the black hole as well as on four 
free functions that measure potential deviations from the 
Kerr metric. This class of metrics includes the Kerr met¬ 
ric as the special case if all deviations vanish. 

My starting point is the Kerr metric which in 
Boyer-Lindquist coordinates is given by the metric ele¬ 
ments 
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In contravariant form, the Kerr metric can be written 
as (e.g., (Hi) 
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General stationary and axisymmetric metrics admit 
two constants of motion, the energy E and axial angular 
momentum and are of Petrov type I. The Kerr met¬ 
ric, however, is of Petrov type D thanks to the existence 
of a third constant of motion, the famous Carter con¬ 
stant which Carter [i^ found by explicitly solving 
the Hamilton-Jacobi equations, 


dSK _l a^pdS^dS^ 
dr Q^fs- U 

As Carter showed, for a Hamilton-Jacobi function of 
the form 

= i/iV - E^t + + S^{r) + Sf{9), (5) 


where r and fj, are the proper time and the rest mass of 
a test particle on a geodesic orbit in this metric, respec¬ 
tively, the Hamilton-Jacobi equations (jd]) are separable in 
all four coordinates making geodesic motion in the Kerr 
metric integrable. I briefly demonstrate the separability 
of the Hamilton-Jacobi equations below, as these steps 


will be essential for my design of a new class of Kerr-like 
metrics. 

From Eq. 0, one obtains the equation 
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E sin^ 9 ^ ^ 

A1 

sV^ry 89 J ■ 


( 6 ) 


After rearranging terms in this equation one can define 
a separation constant such that 



and 



From these expressions, it is obvious that the Hamilton- 
Jacobi equations are separable and one can define the 
Carter constant as 


- (Lf - aE^f. (9) 


In the following, I modify the contravariant Kerr met¬ 
ric in such a manner that the corresponding Hamilton- 
Jacobi equations remain separable. I introduce scalar 
functions /(r), g{9), Ai{r), i = 1,2,5, and Aj{9), j = 
3,4,6, and rewrite Eq. m as 
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where I dehne 


E = E + /(r)+5(0). (11) 


As in the case of the Kerr metric, I can define a 
Hamilton-Jacobi function 

S = - Et + L^(f> + Sr{r) + Se{9), (12) 


where 
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with the corresponding Hamilton-Jacobi equations 


dS _1 dS dS 
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(15) 


Here, 
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(16) 


is the particle’s 4-monientuni. From Eq. (USD 1 obtain 
the equation 
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Substituting the expression m for E and rearranging 
terms 1 can define a separation constant C such that 

(7=^ [-(r^ + a‘^)Ai{r)E + aA 2 {r)L^]'^ 

-p^[r^ + f{r)]-AA,{r)(^^y (18) 

and 

\A3{9)L^-aA4{9)Esm‘^ 9f 

sin 9 

+^l‘^ [a^ cos^ 9 + g{9)] + Aq{9) ■ (19) 

It is now clear that this choice of the functions /, g, 
Hi, i = 1 — 6, preserves the integrability of the Hamilton- 
Jacobi equations for this metric, because all terms which 
depend on either the radius r or the polar angle 9 can be 
separated as in Eqs. (ITSl) and (fTOl) . respectively. Defining 
the Carter-like constant 


I obtain the solutions 
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Sin 6/ 
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Erom here on 1 assume that H 5 (r) > 0 and ^6(0) > 0. 
The former condition is necessary for the regularity of the 
metric as 1 will show in the next section, while asymptotic 
flatness will later require that Ae(9) = 1. 

Setting the partial derivatives of the Hamilton-Jacobi 
function with respect to the constants of motion equal 
to zero, 1 derive the following expressions for the proper 
time and the coordinates t and cj) 
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From the partial derivatives of the Hamilton-Jacobi 
function with respect to the proper time and the coor¬ 
dinates, I obtain the relations between the momenta pa 
and the constants of motion: 
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Finally, using relation (USD, I find the equations of mo- 
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tion of a particle with rest mass fj,: 
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From the contravariant metric specified in Eq. m, I 
obtain the covariant metric 
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This metric reduces smoothly to the Kerr metric in 
Eq. (IT|) if all deviation functions vanish. It is written 
in Boyer-Lindquist-like coordinates, i.e., in spherical-like 
coordinates that reduce to Boyer-Lindquist coordinates 
in the Kerr limit. 

The choice of the deviation functions in Eq. (US may 
not be the most general one, as one might imagine intro¬ 
ducing deviation functions that depend on both the ra¬ 
dius and the polar angle such that the separability of the 
Hamilton-Jacobi equations is still preserved. For mass¬ 
less particles, Eq. m could be multiplied by an arbitrary 
analytic function A^y, 9), which would not spoil the sep¬ 
arability of the Hamilton-Jacobi equations since fi — 0. 
However, it is not obvious whether a similar choice can 
be made for massive particles. Nonetheless, my choice in 
Eq. (nni) is the most general one for deviation functions 
that depend on either the radius or the polar angle. 

Note that the most general stationary, axisymmetric 
metric in general relativity can be written in terms of 
only four functions of the radius and polar angle, if a mi¬ 
nor technical assumption holds. This assumption is satis¬ 
fied in particular for stationary, axisymmetric spacetimes 
that are asymptotically flat and vacuum (see discussion 
in Ref. 0). Thus, one might expect that some of the de¬ 
viation functions /(r) + 9(9), Ai(r) and Aj{9) are either 
trivially related or equal to unity. Indeed, as I will show 


below, requiring that the metric is asymptotically flat re¬ 
duces the number of independent deviation functions to 
four. 

In order to obtain an explicit form of this metric, I 
write the deviation functions Ai{r), i = 1,2,5, as a power 
series in M/r, 
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In Boyer-Lindquist-like coordinates, asymptotically 
flat metrics must be of the form (e.g., [l^, but see dis¬ 
cussion in Ref. (40l| l 
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Expanding the metric given by Eq. (1381) in 1/r, it fol¬ 
lows that oio = 0,20 = 050 = 1 as well as ^ 3 ( 0 ) = 
^ 4 ( 0 ) = Aq{9) — 1. Eurther, I set on = 021 = 051 = 0, 
which defines the parameter M as the mass of the central 
object. Otherwise, the mass would have to be rescaled 
appropriately. This choice also defines the parameter a 
as the spin of the central object. 

The deviation parameters can be further constrained 
in the PPN framework [s^- The general PPN metric can 
be written in the form 


ds^ = —HppN(r)dt^ + Bp-p^{r)dr^ + r^dVt, (43) 
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In general relativity, /IppN = 7ppn = 1- Comparing this 
metric with the expansion of the metric given by Eq. (1381) 
in l/r, I obtain the relations 


2(/3ppN 


7ppn) = £2 — 2ai2 + 


9_i^ 

M2 ’ 


7PPN = 1, 


(47) 

(48) 
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which implies the equation 


^PPN ~ ^ — 2 


£2 — 2a 
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M2 


(49) 


This quantity is tightly constrained by observations [5g : 

|/ 3 ppn - 1 | < 2.3 X ( 50 ) 


In order to avoid any fine-tuning between the parameters 
£2 and ai 2 and the function g(0), I will also set £2 = 
<ai 2 = 5 ( 0 ) = 0. Other choices for these parameters 
that satisfy Eq. 691) could also be made, such as setting 
£2 = 2 q;i 2 and g(0) = 0, but I will exclude them here for 
simplicity. 

Summarizing these results, the newly designed metric 
is given by the elements 

S[A — a‘^A2{r)‘^ sin^ 9] 

[(r2 -I- a2)yli(r) — a’^A2{r) sin^ 0]2 ’ 
a[(r2 a^)Ai{r)A 2 {r) — A]Esin^0 

[(r2-I-a2)Ai(r) — a2A2(r) sin^ 0]2 ’ 

E 

AA 5 (r)’ 

S, 

E sin^ 0 [(r^ -|- a^)'^Ai(r)'^ — a^A sin^ 0] 

[(r2-I-a2)yli(r) — a2A2(r) sin^ 0]2 ^ 


9tt = 

9t(t> — 

9rr — 

9ee = 
9<l><l> ~ 
where 


Ai{r) 

A2{r) 

Mr) 

E 

f { r ) 



(52) 

(53) 

(54) 

(55) 

(56) 


This metric is asymptotically flat, has the correct Newto¬ 
nian limit, and is consistent with the current PPN con¬ 
straints. At lowest order, i.e., truncating the series in 
Eqs. (I52l) -(l56 |) at the first nonvanishing order in the de¬ 
viation parameters, this metric depends on four parame¬ 
ters in addition to the mass M and the spin a: ai 3 , a 22 , 
a52, and £ 3 . 


III. METRIC PROPERTIES 


A. Event horizon, Killing horizon, and ergosphere 


Eirst I find the event horizon. My analysis is similar 
to the one in Ref. [s^. For a stationary, axisymmetric, 
asymptotically flat metric of the kind given by Eq. m, 
the event horizon, if one exists, is located at a radius 

Thor = H{0). (57) 

The function HW) is a solution of the ordinary differen¬ 
tial equation [39|, 

where, from Eq. (l5p. = l/grr and g^^ = l/gee- 

As shown in Ref. , at the poles and in the equatorial 
plane this equation reduces to 

g'^^ = 0 (59) 

due to axi- and reflection symmetry. From the grr and 
ggg elements of the metric in Eq. m, it is obvious that 

E > 0, (60) 

A 5 (r) > 0 (61) 

at all radii r > Phor, because otherwise the metric sig¬ 
nature would change at some radius outside of the event 
horizon which would introduce a singularity at that loca¬ 
tion (see the discussion in the next subsection). There¬ 
fore, Eq. (15^ reduces to the condition 

A = 0 (62) 


both at the poles and in the equatorial plane, which 
means that at these locations the horizon coincides with 
the horizon r+ of a Kerr black hole. 


r += HK = M + v'M 2 - a 2 . ( 63 ) 

As in Ref. [s^, here I use both notations r+ and Hk for 
the Kerr horizon. While the former is typically used, the 
latter is often used in the numerical relativity literature. 

For arbitrary angles 9, the horizon equation (1551) takes 
the form 


1 

E 


AA^iH) + 



= 0 , 


which reduces to the equation 


(64) 


(65) 


In this section, I analyze some of the properties of the 
newly designed metric and calculate the location of the 
event horizon. I likewise determine conditions which the 
deviation functions in Eqs. 61I)”6ni) must satisfy so that 
the exterior domain is regular. 


because E > 0. Since also A^{H) > 0, it is obvious that 
the Kerr horizon given by Eq. dM]) is always a solution 
of Eq. dM]). Furthermore, should a different solution to 
the horizon equation exist, it must lie inside of the ra¬ 
dius r+ at angles 0 < 9 < 7 r/ 2 , 7 r /2 < 0 < tt, because 













the second-term in Eq. (l65|) is non-negative and ^5 > 0, 
which implies that a solution for which A ^ 0 must have 
A < 0 and, therefore, r < r+. From this analysis I con¬ 
clude that the event horizon of the metric in Eq. (EU is 
the Kerr event horizon. 

The Killing horizon of the metric given by Eq. 0 is 
located at the roots of the equation 


9tcf) 9tt94>4> 


AE2 sin^ 0 


[(r2 -I- a^) Ai{r) — a'^A2{r) sin^ 0] 


r, (66) 


which is equivalent to the equation 


A = 0 (67) 

as long as the denominator of Eq. (l66l) is nonzero at and 
outside of the Killing horizon. I will return to this re¬ 
quirement below. Consequently, the Killing horizon co¬ 
incides with the event horizon. This suggests that a gen¬ 
eralized version of Hawking’s rigidity theorem Q may 
hold for this metric. 

The ergosphere is determined by solving the equation 

9tt = 0, (68) 


which reduces to the equation 

A-a^A2{r)'^sm^0 = O. (69) 


The shape of the ergosphere, therefore, depends on the 
spin and on the function A 2 (r). 

At lowest nonvanishing order, Eq. ((M)) reduces to the 
equation 

A - {r^ + a22M‘^f 8111^0 = 0, (70) 
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FIG. 1. Event horizon and ergosphere in the a; 2 -plane of a 
black hole with spin |a| = 0.9M for several valnes of the pa¬ 
rameter 022 . The event horizon coincides with the horizon of 
a Kerr black hole and depends only on the spin a. For increas¬ 
ing values of the parameter 022 , the extent of the ergosphere 
increases at angles 0 < 0 < vr. At the poles, the event hori¬ 
zon and ergosphere coincide irrespectively of the value of the 
parameter Q 22 . 


For the determinant to be negative definite, I impose the 
requirement 


which has to be solved numerically. In Fig. [1] I plot the 
event horizon and the ergosphere of a black hole with a 
spin |o| = 0.9M for several values of the parameter a ;22 
in the a; 2 :-plane, where x = + a? sin 0 , 2 ; = rcos 0 . 

For increasing values of the parameter a 22 , the extent of 
the ergosphere increases. 


B. Regularity of the exterior domain 


Following the analysis in Ref. [s^, I analyze the exte¬ 
rior domain of the metric given by Eq. m for the exis¬ 
tence of any singularities or pathological regions, which 
might have been introduced by the deviation functions 
in Eqs. (I5^ -(I55 )) . Arbitrary deviations from the Kerr 
metric can lead to a violation of Lorentzian signature or 
the existence of closed timelike curves. In the former 
case, the determinant of the metric is no longer negative 
definite and, in the latter case, the {(p, (p) element of the 
metric is negative. 

The determinant of the metric in Eq. m is given by 
the expression 


det (gap) = 


sin^ 0 

A5(r) [Ai(r)(r2 -|- a^) — a?A2(r) sin^ 0]^ 

(71) 


Ai(r){r^ + a^) - A2(r) sin^ 0^0, (72) 


in addition to the conditions in Eqs. dm]) and dna, which 
must be fulfilled everywhere on and outside of the event 
horizon. 

At the lowest order in the deviation parameters, I can 
write these conditions in the following manner: 


a^r(r^ -|- 0:22Af^) sin^ 0 — r^(r^ -|- a?) 
^ M3(r2 + a2) 


<^52 > — 


£3 > - 


(M +VM2-a2)" 
{M + yf - a'^f 


(73) 

(74) 

(75) 


In these expressions I replaced the radius with the Kerr 
horizon radius and, thereby, obtain a lower bound 
on the deviation parameters, which is valid at all radii 
r > r+. 

In order to exclude the existence of closed timelike 
curves in the exterior domain, I derive from the (</>, p) 
element of the metric the additional requirement 


Ai(r)2 -I- - a^AsA? 0 > 0 


(76) 
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or, at lowest order, 


ai3 > — 


{M + y/M^-o?f 


(77) 


on and outside of the horizon, where in the last equation 
I, again, inserted the Kerr horizon radius r = r+. Com¬ 
bining Eq. ([77| with the condition given by Eq. ([731). I 
obtain the requirement 


0122 > ~ 


(M-b VM2-a2)^ 


(78) 


Therefore, Eqs. dill), dZSl), dZZl), and (l78ll define the al¬ 
lowed ranges of the deviation parameters in the lowest- 
order metric as a function of the spin. I plot these func¬ 
tions in Fig. [2j At higher values of the spin, the lower 
limit on the deviation functions is more tight than at 
smaller spin values. Note, however, that at maximal spin 
a = ±M the lower limit on all four parameters is — 1, 
which means that even in this case both positive and 
negative deviations from the Kerr metric can be studied. 
This is different in, e.g., the metric designed in Ref. [s^, 
where the metric properties can change drastically for 
positive deviations larger than a lower bound that ap¬ 
proaches zero in the limit a —>■ ±M [ssl 

I verified numerically the absence of curvature singu¬ 
larities at the horizon located at r = r+ for a large 
number of different values of the spin in the range 
— 1 < a/M < 1 and of the deviation parameters 
£ 3 , ai 3 , a 22 , a 52 in the allowed parameter space by eval¬ 
uating the Kretschmann scalar. 


K = (79) 

at polar angles 0 = 7r/2 — 57rn/21, n = 0,1,2, where 
R°‘pryS Riemann tensor. For this analysis, I chose 

equidistant values of the spin in steps of Aa = O.IM and 
values of the deviations parameters ranging from —3 to 5 
in steps of 1 neglecting any parameter combinations that 
lie in the excluded part of the parameter space. In all 
cases, the event horizon is the location of a coordinate 
singularity as in the Kerr metric. In Sec. El I will con¬ 
struct an explicit coordinate transformation that removes 
the coordinate singularity at the event horizon, proving 
that it is indeed the location of a coordinate singularity. 
I, therefore, conclude that the metric in Eq. m harbors 
a black hole with a regular exterior domain. 


IV. KEPLERIAN FREQUENCY, ENERGY, 
ANGULAR MOMENTUM 


In this section, I derive expressions for the dynamical 
frequencies, energy, and axial angular momentum of a 
particle on a circular equatorial orbit. My derivation 
follows closely the one in Refs. Then I proceed to 

calculate the radius of the ISCO. 



FIG. 2 . Lower bounds on the deviation parameters aia, 022, 
052, and £3 of the metric in Eq. at the lowest nonvan¬ 
ishing order of the deviation parameters. For values of the 
deviation parameters greater than these bounds, the event 
horizon and the exterior domain are regular, i.e., free of cur¬ 
vature singularities and closed timelike curves. 


The geodesic equations, 

dx^ dx^ 

dr'^ dr dr ’ 


(80) 


where are the Christoffel symbols, can also be written 
in the form 

d / \ 1 dgp-y dx^ dx'^ 

dr dr J 2 dx^ dr dr 

Due to axi- and reflection symmetry, dr/dr = dd/dr = 
d^r/dr^ = 0 for particles on circular equatorial orbits. 
Therefore, the equation that governs particle motion in 
the radial direction reduces to the relation 


drgtd'^ + 2drgt<j>i(l) + 9^500^^ = 0, (82) 


where dr = d/dr and the dot denotes the derivative d/dr. 
From this equation, the Keplerian frequency = (p/t = 
P 4 >/pt as observed at radial infinity can be expressed in 
the form 


„ drPtrp i \/ {drgttj))^ /'oo^ 

^-) (o3j 


where the upper sign refers to prograde orbits, while the 
lower sign refers to retrograde orbits. Comparing this 
expression with the metric in Eq. dsn, it is obvious that 
the Keplerian frequency is independent of the function 
A^{r), because it is independent of the metric element 

Qrr ■ 
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FIG. 3 . Dependence of the Keplerian frequency = c^D,^/2nGM on the deviation parameters es, aia, and 022- Each panel 
shows the Keplerian frequency v,j, as a function of radius for a black hole with mass M = IOMq and spin a = 0 . 8 M for different 
values of one deviation parameter while setting the other two equal to zero. At a given radius, the Keplerian frequency increases 
for decreasing values of the parameters es and 022 and for increasing values of the parameter a 13. The dot denotes the location 
of the ISCO. 


Eq. (IMll assumes that drg^^ ^ 0. At least in the case 
of the lowest-order metric, for large negative values of the 
deviation parameter ai 3 near the lower bound given by 
Eq. (l77l) . the term can be either zero or negative 

near but outside of the event horizon; otherwise it is al¬ 
ways positive. If the term drg^,^, vanishes, the Keplerian 
frequency is simply = —drgu/‘^drgt<t> from Eq. (15^ . 

In Eig. |3l I illustrate the dependence of the Keplerian 
frequency on the three lowest-order deviation parameters 
£ 3 , ai 3 , and a 22 . I plot the Keplerian frequency for a 
black hole with mass M = IOMq and spin a = 0.8M 
as a function of the radius for different values of one 
of these parameters while setting the other two equal to 
zero. At a given radius, the Keplerian frequency increases 
for decreasing values of the parameters £3 and q ;22 , while 
the Keplerian frequency increases for increasing values of 
the parameter 013 . 

From the equation 




(84) 


evaluated in the equatorial plane, where is given by 
Eq. (fT6)) . I obtain the effective potential 

hefF(?’) = 2^* (jirr'f' + gOB^ ^ 

= - 2g*^EL, + g^^Ll + g^). (85) 

In this expression, I substituted the constants of motion 
E and given by Eqs. O-dSH) for the momentum 
components pt and Circular equatorial orbits are 
governed by the equations 


Eeff(r) = 0 , 

dVesir) 


( 86 ) 

(87) 


In terms of the constants of motion, I write the Kep¬ 
lerian frequency as 


^ _ gt4,E + guLz 

Pt g4><pE -f gt(f/Ez 


( 88 ) 


Combining this equation with Eqs. (|M1) and (IMll . I can 
solve for the energy and axial angular momentum and 
obtain the expressions 


E 


gtt + gttji^tj) 


^ ^ 5tt ‘2gt4>^4> 500^^ 

gtip 4 “ grf>rf>^4' 




^ sj 5tt ‘^gt<f>^<f> 300^^ 


(89) 

(90) 


where, again, the upper and lower signs refer to prograde 
and retrograde orbits, respectively. These expressions are 
likewise independent of the deviation function A^{r). 

In the metric designed in Ref. [s^, for extreme val¬ 
ues of the deviation parameters and at certain radii very 
close to the ISCO, sign changes can occur in the expres¬ 
sions of the energy and axial angular momentum so that 
these quantities are smooth at these radii; see Ref. [ 2 ^ 
for a detailed discussion. In principle, a similar behavior 
in these expressions of the metric designed in this paper 
could arise, if the discriminant of the square root in the 
denominator of Eqs. (1591) and o has roots. Empiri¬ 
cally, by evaluating these functions for different values of 
the deviation parameters, I find that such roots exist at 
least in the case of the lowest-order metric. However, the 
energy and axial angular momentum are still smooth at 
these roots and I have seen no evidence for discontinu¬ 
ities of the first derivatives of the energy or axial angular 
momentum outside of the ISCO. Should these nonethe¬ 
less exist, then the opposite sign of the square roots in 


dr 
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FIG. 4 . Dependence of the energy E of a, particle with rest mass ^ on a circular equatorial orbit around a black hole with spin 
a = 0.8M as observed at infinity on the deviation parameters £3, 013, and 022- Each panel shows the energy as a function of 
radius for different values of one deviation parameter while setting the other two equal to zero. At a given radius, the energy 
increases for decreasing values of the parameters £3 and 0:22 and for increasing values of the parameter 013. The ISCO is located 
at the minimum of the energy and is denoted by a dot. 





Radius (M) Radius (M) 


FIG. 5 . Dependence of the axial angular momentum Lz about the black-hole spin axis of a particle with rest mass /r on a 
circular equatorial orbit around a black hole with spin a — 0 . 8 M as observed at infinity on the deviation parameters £3, 013, 
and a22. Each panel shows the axial angular momentum as a function of radius for different values of one deviation parameter 
while setting the other two equal to zero. At a given radius, the axial angular momentum increases for decreasing values of the 
parameters £3 and 022 and for increasing values of the parameter 0:13. The dot denotes the location of the ISCO. 


Eqs. (l89l) and (1^ has to be chosen at radii smaller than 
the root of the discriminant. 

In Figs, m and [5l I plot the energy and axial angular 
momentum versus the radius r for different values of the 
lowest-order deviation parameters. As before, on varying 
one of these parameters, I set the other two equal to 
zero. At a given radius, both the energy and the axial 
angular momentum increase for decreasing values of the 
parameters £3 and 022 and for increasing values of the 
parameter ai 3 . 

I calculate the location of the ISCO from the equation 


In Fig. | 6 l I plot contours of constant ISCO radius versus 
the spin for different values of the deviation parameters. 
Again, on varying one of these parameters, I set the other 
two equal to zero. 

The dependence of the ISCO on these parameters is 
more complex. As shown in Fig. [BJ at a fixed value of 
the spin, the location of the ISCO decreases for increasing 
values of the parameter £3 for values of the spin a ^ 
0.8M. For values of the spin a > 0.8M, the ISCO radius 
increases for increasing values of the parameter £3. At a 
spin a ~ 0.8M, the ISCO is practically independent of 
the parameter £3. The location of the ISCO increases for 
increasing values of the parameter ai 3 except for values of 
the spin a ^ 0.7M and values of the parameter ai^ which 
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FIG. 6. Dependence of the ISCO radius on the deviation parameters ea, Q13, and 022- Each panel shows contours of constant 
ISCO radius as a function of the spin and one deviation parameter while setting the other two equal to zero. At a fixed value 
of the spin, the location of the ISCO can either increase or decrease for increasing values of each deviation parameter and, in 
some cases, be practically independent of each deviation parameter. In the green shaded region of the central panel, the energy 
has two local minima and the ISCO is located at the outer radius where these minima occur. In the red shaded region of the 
central panel, circular equatorial orbits do not exist at radii r 2 . 5 M and the ISCO is located at the outer boundary of this 
radial interval. The black shaded region marks the excluded part of the parameter space. 


are very close to the boundary as defined in Eq. G?]), 
where the ISCO becomes practically independent of the 
parameter 0 : 13 . If the ISCO is located at a radius nsco ^ 
3M, it depends only weakly on the parameter ai^. 

For values of the parameter ai 3 ^ —5.7, the energy E 
can have two local minima. In this case, stable circular 
orbits exist only at and outside of the radius where the 
energy has a local maximum between its two local min¬ 
ima, because the radial epicyclic frequency is imaginary 
otherwise (see the discussion below), and the ISCO is lo¬ 
cated at the outer radius where the minima occur. For 
values of the parameter ai 3 ^ — 6 . 2 , both the energy and 
Keplerian frequency become imaginary in a small range 
of radii r ^ 2.5M and the ISCO simply lies at the outer 
boundary of this radial interval even if the energy does 
not have a local minimum at this location. In certain 
cases, the Keplerian frequency vanishes approximately 
at the ISCO. I illustrate this behavior of the energy and 
Keplerian frequency in Fig. [71 

Finally, the ISCO radius is independent of the param¬ 
eter 022 if a = 0 and depends only weakly on the pa¬ 
rameter 022 if |a| ~ 0 or if nsco ^ 3M. If a > 0, the 
location of the ISCO decreases for increasing values of the 
parameter 022 , while the location of the ISCO increases 
for increasing values of the parameter 022 if a < 0 . 

In order to derive expressions for the radial and vertical 
epicyclic frequencies fir and Og, I write Eq. (I5S1) in the 
forms 


1 

2 

1 

2 



Veff 

5rr(p*)^ 

I4ff 

geeip^Y 



(92) 

= Kff- 

(93) 


Now I introduce small perturbations 5r and 59 and 
take the derivative of Eqs. ([Ml) and (IMI) with respect to 
the coordinate time, which yields the equations 




_ cPV^ 
dE 

d^{S9) 
dt"^ 


eff 


5r, 




(94) 

(95) 


From these expressions, I derive the radial and vertical 
epicyclic frequencies as 


0 " = — 


0 ? = — 


|2wr 

(96) 

'Xff 

d02 > 

(97) 


where the second derivatives are evaluated at r = tq. 
These expressions are lengthy and I do not write them 
here explicitly. Since the metric element grr occurs only 
in Eq. (|92ll and not in Eq. (IM|l . the vertical epicyclic 
frequency is independent of the deviation function 
while the radial epicyclic frequency depends on all four 
deviation functions. 

In Figs. |5| and ini I plot the radial and epicyclic frequen¬ 
cies for a particle on a circular equatorial orbit around 
a black hole with mass M = IOMq and spin a = 0.8M 
for different values of the deviation parameters. At a 
given radius, both the radial and the vertical epicyclic 
frequency increase for increasing values of the parame¬ 
ters £3 and Q !22 and for decreasing values of the parame¬ 
ter ai 3 . The radial epicyclic frequency also increases for 
increasing values of the parameter 052 . 

For all values of the deviation parameters in the al¬ 
lowed part of the parameter space as shown in Fig. [7] up 
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FIG. 7 . Energy (top) and Keplerian frequency (bottom) of 
a particle on a circular equatorial orbit around a black hole 
with spin a = 0 . 4 M for several values of the parameter Q13; 
all other deviation parameters are set to zero. For large nega¬ 
tive values of the parameter 0:13, the energy has a second local 
minimum inside of the ISCO radius (labeled by a dot) with 
a higher value of the energy than at the ISCO. The Keple¬ 
rian frequency decreases at radii just outside of the ISCO and 
vanishes approximately at the ISCO radius. For even more 
negative values, both the energy and the Keplerian frequency 
are imaginary at radii r ~ 2 . 5 M and circular equatorial orbits 
no longer exist. In the latter case, the ISCO lies at the outer 
boundary of this radial interval. 


to values of at least - 1 - 10 , the radial epicyclic frequency 
always vanishes at some radius outside of the event hori¬ 
zon, which coincides with the ISCO unless the energy E 
has more than one local minimum (c.f., the green and 
red shaded regions in Fig. 15]). Circular equatorial orbits 


are radially unstable inside of this radius and plunge into 
the black hole. For values of parameter 013 in the green 
shaded region of the central panel in Fig. [SJ the radial 
epicyclic frequency vanishes at the radius where the en¬ 
ergy has a local maximum. In this case, the ISCO is 
located outside of this radius and the vertical epicyclic 
frequency has a minimum near the ISCO radius. I il¬ 
lustrate this property of the radial and vertical epicyclic 
frequencies in Fig. [TOl where I plot these frequencies for 
a black hole with mass M = IOMq and spin a = 0.4M 
for a value of the parameter 033 = — 6 . 2 . 


V. KERR-SCHILD FORM 

In this section, I derive expressions for the principal 
null congruences of the metric given in Eq. (EU and con¬ 
struct a transformation of the Boyer-Lindquist-like co¬ 
ordinates to Kerr-Schild-like coordinates. My derivation 
is similar to the corresponding calculations for the Kerr 
metric (c.f., e.g., Ref. (b^). 

For null geodesics, it is convenient to introduce the 
parameters 

(98) 

‘^=§ 2 ’ ( 99 ) 

so that the functions R{r) and 0(0) in Eqs. (|23)) and 

(1^51) can be written in the form (note that /i = 0 and 
AsiB) = Aii0) = Ae{B) = 1) 

= [(r^ -b a^)Ai{r) - 0^2(0^]^ 

-A[(e-a)2 + ry], (100) 

=r]+{^-af - ^^-i^-asin^Bf. ( 101 ) 

-c/ Sin 0 

Since from Eq. 0(^) > 0, 

0+(e-a)'>O, (102) 

where equality holds if and only if 

B = Bq = const., (103) 

^ = asin2 0. (104) 

In this case, 

77 = — cos"^ 00- (105) 

These expressions for the parameters ^ and 77 are identical 
to the ones for the Kerr metric. 

Using the equations of motion, Eqs. (El-dSZl), and 
Eqs. (|104l) and (11051) . I obtain the photon equations of 
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FIG. 8 . Dependence of the radial epicyclic frequency Vr = c?Q.r/2 -rGM on the deviation parameters £ 3 , 013 , 022 , and 052 . Each 
panel shows the radial epicyclic frequency as a function of radius for a black hole with mass M = IOMq and spin a = 0.8M 
for different values of one deviation parameter while setting the other two equal to zero. At a given radius, the radial epicyclic 
frequency increases for increasing values of the parameters £ 3 , 022 , and 052 and for decreasing values of the parameter 013 . 
The ISCO is located at the radius where the radial epicyclic frequency vanishes. 


motion 


tions are given by the vectors 


dt 

dA 


dr 

dA 


de 

d(l) 



a'^)Ai{r) 

A 


(r^ + aF‘)Ai{r) — a‘^A2{r) sin^ Oq 
S 

±rE = ±^A5ir) 

(r^ + a^)Ai{r) — a^A2{r) sin^ Oq 


0 , 

A 

(r^ + a^)Ai(r) -a^A2(r)sin^ Op ^ 


where A is an affine parameter. 

Accordingly, setting E = 1, the principal null direc- 


11 = {l\±l\0,l’^) . (107) 

In order to remove the coordinate singularity of the 
metric in Eq. m written in Boyer-Lindquist-like coor¬ 
dinates located at the event horizon r_|_, I perform a trans¬ 
formation to Kerr-Schild-like coordinates, which can be 
defined using either the outgoing {+E) or ingoing {—F) 
principal null direction in Eq. (11071) . Defining the trans¬ 
formations 


dtKS = dtBi, T -j^dr^L, 

drKS = EydrBL, 
dOKS = d0Bh, 

F 

d(l)KS = dipBhT-jydvBh, ( 108 ) 

the transformed principal null vectors (setting E = 1) 
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FIG. 9 . Dependence of the vertical epicyclic frequency ve = c^Q,e/ 2 'rGM on the deviation parameters £3, ais, and 022- Each 
panel shows the vertical epicyclic frequency Vr as a function of radius for a black hole with mass M = lOM© and spin a = 0 . 8 M 
for different values of one deviation parameter while setting the other two equal to zero. At a given radius, the vertical epicyclic 
frequency increases for increasing values of the parameters es and 0:22 and for decreasing values of the parameter 013. 



FIG. 10 . Radial and vertical epicyclic frequencies of a particle 
on a circular equatorial orbit around a black hole with mass 
M = lOM© and spin a = 0 . 4 M for a value of the parameter 
ai3 = —6.2; all other deviation parameters are set to zero. 
The black dots marks the location of the ISCO. The radial 
epicyclic frequency vanishes at the radius that corresponds to 
the local maximum of the energy E (c.f., top panel in Fig. [T]). 
The vertical epicyclic frequency has a local minimum at the 
ISCO. 


are given by the expression 

;j = (o, ±1,0,0), (109) 

which is the form of the principal null directions of the 
Kerr metric in standard Kerr-Schild coordinates. Here, 
the upper/lower sign refers to the transformation in 


Eq. (11081) with the upper/lower signs and the subscripts 
“KS” and “BL” stand for “Kerr-Schild-like” and “Boyer- 
Lindquist-like”, respectively. Note, however, that the 
transformation of the radius makes this transformation 
rather cumbersome in practice and it may not be possible 
to perform it explicitly. 


For this reason, I define an alternative transformation 
to Kerr-Schild-like coordinates, which drops the transfor¬ 
mation of the radius in the transformation in Eq. (11081) 
and which is given by the relations 


dtKS 

dvKS 

ddKS 

d(j)KS 


dtBL 


drBh, 

ddBL, 

dcjiBh 


2MrAi (r) 

^\/A5(r) 


drBL, 


aA2(r) 


drBL- 


( 110 ) 


In these expressions, I used the definitions of the null 
vector in Eq. (11071) given by Eq. (11061) with a slight mod¬ 
ification of the component I*, where I replaced the factor 
(r^±a^)/A by 2Mr/A as in McKinney and Gammie (d^ . 
A transformation to Kerr-Schild-like coordinates can be 
defined for both of these choice. 


The metric in the above Kerr-Schild-like coordinates is 
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then given by the elements 

pjg 5i[A — a^A2(r)^ sin^ 0] 

9tt = 


KS _ 

i/tr “ 


-a^A^ir) sin^0}, 


F 


{Ai{r)[2Mr + a^A 2 {ry sin^ 9] 


KS _ oS[(r^ + ayAi{r)A 2 {r) - A] sin^ 9 
9t<j> — 


KS ^ 

A 5 {r)F 

— 2 a‘^A 2 {r) sin^ 0 }, 

KS qS sin 9 

9r<p 


F 

{Ai(r)[A + 4Mr + a?A 2 {ry sin^ 9\ 

'■e}. 

[(r^+c?)AyrfA 2 {r) 


^/^F 

+2MrAi{r) — a^A 2 {r) sin^ 0], 


9fe^ = S, 


KS 


E [(r^ + Ai (r)^ — A sin^ 9'\ sin^ 9 


F 


, ( 111 ) 


where 


F= [{r^ + ayAi{r) - a^A 2 {r)sm^ 9]^ . (112) 


In the Kerr case, Ai{r) = A 2 {r) = A^r) = 1, /(r) = 0, 
and this metric reduces to the Kerr-Schild metric in the 
form given in Eq. (4) of McKinney and Gammie [i^ . 

From the metric elements in Eq. (11111) it is clear that 
the metric in Kerr-Schild-like coordinates no longer has a 
coordinate singularity at the radius given by Eq. (1551) . 
Evaluating the (r, r) and {9,9) elements of the contravari- 
ant metric and inserting them into the horizon equation 
(l58)) . I obtain an equation that is identical to Eq. (1651) . 
Therefore, the metric in Kerr-Schild-like coordinates still 
harbors an event horizon at the radius r+. 

The Kerr-Schild-like form of the metric in Eq. (11111) 
is important for fully relativistic magnetohydrodynamic 
simulations of accretion flows in this metric, because 
these typically require the presence of an event horizon 
without any coordinate singularities at the horizon or in 
its vicinity both inside or outside of the horizon (see, e.g. 
Ref. [1^). Otherwise, matter in the simulation might get 
trapped at the horizon leading to an unphysical matter 
accumulation. In particular, the lapse 


N = 


1 

V-9ks 


(113) 


has to be positive across the event horizon. 

The lapse of the lowest-order metric in Kerr-Schild-like 
form is given by the expression 


N = 


^5(j,3 _|_ ggJVfS _|_ Q'j 


y (r^ -I- 2Mr -|- a^)(r3 -|- sin^ 9 ’ 

(114) 

which is indeed positive across the horizon in the allowed 
part of the parameter space [recall Eq. (1551) ]. In Fig. ITTl 



Radius (M) 



Radius (M) 


FIG. 11 . Lapse versus radius for a black hole with mass M 
and spin a = 0 . 8 M in the equatorial plane for different values 
of the parameters (top) ea and (bottom) Q13. In each panel, 
only one parameter is varied, while the other one is set to 
zero. The lapse is positive across the radius that marks the 
location of the event horizon (vertical dashed line). 


I plot the lapse for a black hole with mass M and spin 
a = 0.8M in the equatorial plane for different values of 
the parameters €3 and ais. 

Accretion flows in three-dimensional fully relativistic 
magnetohydrodynamic simulations are greatly affected 
by the amount of frame-dragging near the black hole, 
while the location of the ISCO is only of marginal impor¬ 
tance. In order to assess the amount of frame-dragging of 
the lowest-order metric, I expand the (t, (f>) component of 
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the metric in Eq. (Illll) in 1 /r and obtain the expression 


9t(p 


2aM sin^ 6 022 -^ 

- 1 + ^^ 

r I 2r 

2a? cos^ 0 — aisM'^ 



(115) 


match metrics of this form with the metric in Eq. m, 
which I expand to linear order in the deviation parame¬ 
ters. 


A. The modified gravity bumpy Kerr metric 


The leading nonvanishing order of the deviations from 
the Kerr metric is of order 1/r^, while the next-order 
correction of the Kerr part is only of order 1/r^. Con¬ 
sequently, frame-dragging can be strongly affected by 
the presence of non-Kerr deviations, and positive values 
of the parameters 022 and 013 enhance the amount of 
frame-dragging. I illustrate the dependence of the (f, (j)) 
element of the metric in Kerr-Schild-like form on the pa¬ 
rameters ai 3 and 022 in Fig- HH 

VI. MAPPING TO OTHER METRICS 

In this section, I derive the explicit mapping of the 
metric given in Eq. m to other known metrics from 
the literature. In particular, these include known four¬ 
dimensional, analytic black hole solutions in modified 
theories of gravity. These metrics can be written as 

T (116) 

where is the modification of the Kerr part of the 
metric, which depends on one or more parameters. I will 


The metric designed in Ref. 0 (the “modified gravity 
bumpy Kerr metric”, hereafter labeled “MGBK”) derives 
from the most general stationary, axisymmetric metric in 
Lewis-Papapetrou form with the additional requirement 
that it possesses three independent constants of motion 
for small deviations away from the Kerr metric, where 
the third, Carter-like constant is quadratic in the mo¬ 
mentum. This metric is defined in terms of a coupled 
set of integro-differential equations, which (in the de¬ 
formed Kerr parametrization of Ref. 0 ) are listed in 
Appendix A. Note that this metric depends on four de¬ 
viation functions as is the case of the metric in Eq. m, 
which are denoted 'yi{r), i = 1,3,4, and © 3 ( 0 ). 

Making an ansatz for the (t, (f>) component of this met¬ 
ric in the form of a power series in l/-\/S, 

/ A/fn \ 

(117) 

Gair and Yunes 0 solved the set of integro-differential 
equations and obtained the remaining metric elements: 


t MGBK 

'ht 


^MGBK 


tMGBK 


t MGBK 


_ ^ ^2 7 MGBK 

M a 


a (r'^ + a?) ji?Is. 2a^r^Asin^6* 


Asin^ 


2M 


7i 


Pi dr Pi 

EA^ sin^ 6* dyi a (S + 2a^Mr sir? 9) sir? 6 


-71 


M 


-73 


A^(E — 4Mr) sin^ 6* ^74 


0 , 


2M F 
S 71 
' A ’ 

(r^ -I- a^)^ 


dr 2M 


a 


dr 2M 


a 


dr 


tMGBK 


-2^71- 


2(r^ -|- a?) 


T MGBK 

^tcf> 


2A^ sin^ 9 


-73 


2A2 

^- 74 , 


2AP, 


E a 


7r74 


(118) 


where 


-2 _ 2 2 2/1 

1 = r — a cos 9, 


©3(0) = 0, 



A =1,4, 


^ 00 


73 = - 


. 13 ,r. 


n—0 



(119) 

( 120 ) 

( 121 ) 

( 122 ) 


The functions i = 1 — 4, are polynomials in r that can 
be found in Appendix A of Ref. [ij and the coefficients 
ht(ji,n are given in Ref. [l5j |. Additional simplifications [lB | 


allow one to set 

71.0 = 71.1 = 73,0 = 73,2 = 74,0 = 74,1 = 0. (123) 

Gair and Yunes 0] gave an explicit version of this 
metric expanding all metric components in power series 
in 1/r, 


= (124) 

n ^ ^ 


where the coefficients are given in Eqs. (27)-(30) 

of Ref. 01 . Note that while current PPN constraints 
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FIG. 12 . Element of the lowest-order metric to 0 {l/r^) 
versus radius for a black hole with mass M and spin a = 0 . 8 M 
in the equatorial plane for different values of the parameters 
(top) ai3 and (bottom) 0:22 ■ In each panel, only one param¬ 
eter is varied, while the other one is set to zero. This ele¬ 
ment indicates the amonnt of frame-dragging, which depends 
strongly on the deviation parameters in the vicinity of the 
event horizon. The vertical dashed line marks the location of 
the event horizon. 


were taken into account in Ref. 113 up to 0 {l/r), the 
htt,2 coefficient in Eq. ( 27 ) of Ref. [l^ is likewise tightly 
constrained by current PPN experiments (c.f., Eq. (| 4 ^ 
and Ref. [s^l) and 


71,2 = 73,1 = 74.2 = 0 ( 125 ) 

should be chosen so that htt,2 = 0 . In principle, the 
choice 71,2 = —274,2, 73,1 = 0 is sufficient so that the 
htt,2 coefficient complies with the PPN constraints. This 


choice, however, would introduce an undesirable finetun- 
ing between the parameters 71,2 and 74 , 2 , which was also 
avoided in some of the other simplifications in Ref. [l^ 
as listed in Eq. (11231) . 

With this minor adjustment, I obtain from the metric 
elements and the mapping 

asn = 7i,n, n>2, (126) 

e„ = 0, n > 3, (127) 

as well as from the elements the mapping 

2“- t) = 

n—3 ^ ^ ' 

-b -b oV(3r-b 4M) 

-ba2(A - 2Mr) cos 26»]/i“](^28) 
°° / M \" 1 

n^2 ^ '' 

+2{i:-2Mr)csc^0h^^^^], (129) 

where the right-hand side of Eqs. (I128I) " (I129I) has to be 
expanded in 1/r. Should any terms in this expansion de¬ 
pend on the polar angle 6 , the coefficients of these terms 
have to be set to zero by an appropriate choice of the de¬ 
viation parameters 7fe,n, A: = 1,3,4, n > 2, because the 
left-hand side in these equations are functions of radius 
only. As I will show below, no such angular terms oc¬ 
cur at least up to order n = 5. Note that the remaining 
matching of the metric element with the devia¬ 

tion from the (c^, (f>) element of the metric in Eq. (1511) 
linearized in the deviation parameters serves here as a 
consistency check of the mapping of the deviation param¬ 
eters in Eqs. (11281) ^ (11291) . which may impose additional 
requirements on certain of the parameters 'yk,n- This is 
likewise not the case at least up to 0 {l/r^). 

The coefficients are given in Ref. [l^ up to 

0{\/r^), and up to this order I find for the coefficients 
of Eqs. (|128l) - p29l) the explicit relations 

= 74 , 3 , 

= - 274,3 + 74 , 4 , 

= —274,4 + 74,5) 

M 

= -73,3, 

a 

M 

= -( 273,3 - 73 , 4 ), 

a 

a \ 

= 7773,3 -( 273,4 - 73 , 5 ), 

M a 

a M 

= (^73,4 -—( 273,5 - 73,6). (130) 

These equations suggest a general mapping of the form 

C^5n _ r, I 

CTln 2 — ^74,n—1 “b 74,n, 

a2n -^ - - 73,n-|-l), (131) 
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where the first equation holds for n > 3, while the second 
equation holds for n > 2 [note Eqs. (11231) and (I125p ]. I 
have not investigated this mapping at orders n > 5. If 
this mapping is valid at all orders n, then, together with 
Eqs. (11261) and (11271) . the linearized form of the metric in 
Eq. ((5T|| and the metric of Ref. [l^ in the above form [l^ 
are equivalent. Eq. (11271) reduces the number of deviation 
functions of the metric in Eq. m for this mapping to 
three, which is in accordance with the choice of Gair and 
Yunes [l^ to set 03 ( 0 ) = 0. 

Due to the implicit form of some of the elements in the 
general metric of Ref. 0 in Eq. m I do not attempt 
to map it to the linearized form of the metric in Eq. m 
or in Eq. ([38|) . Should it turn out that these two met¬ 
rics cannot be mapped exactly, then it must be possible 
to further generalize the metric designed in this paper, 
because the metric in Eq. dSB) was obtained from the 
most general stationary, axisymmetric metric in Lewis- 
Papapetrou form which admits three constants of motion 
for small deviations from the Kerr metric and because the 
Carter-like constant in both metrics is quadratic in the 
momentum (c.f., Eq. (1321) in this paper and Eq. (38) in 

Ref. 0). 


B. Einstein-Dilaton-Gauss-Bonnet Gravity 


only the (t, (f>) component of the metric is modified, which 
is given by the expression 


- CS 
H<j> 


a 



12M 

7r 


27M^\ 

10r2 ) 


(135) 


In this case, the mapping is 


^24 = gCcS, 

a25 = Y^Ccs, 

27 

a 26 = TxCcs- 

16 


(136) 

(137) 

(138) 


All other deviation parameters vanish. Note that the 
metric to O(a^) found in Ref. 0 is not integrable and, 
thus, cannot be mapped to the metric in Eq. m- 


D. Braneworld black holes 

One class of metrics that cannot be related to the met¬ 
ric in Eq. dSlD via a simple mapping is the rotating black 
hole solution in Randall-Sundrum-type braneworld grav¬ 
ity 01 , which was found in Ref. 0| . This metric is given 
by the elements 0| 


Static black holes in gravity theories described by La- 
grangians modified from the standard Einstein-Hilbert 
form by scalar fields coupled to quadratic curvature in¬ 
variants were investigated in Ref. 0]. In these solutions, 
the relevant component of the metric deformation is given 

by 


7 EDGE _ 0^3 

rr ^^2^2Jg(^)2 

/ M 52 2M3 

1 H-f - T- H-5“ 


5 r4 3 r5 ’ 


where is the coupling constant of this theory and 
K = l/(167r). 

The mapping is, then, given by the equation 


“S’! 
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mV 




a^M'^ 


15Kr®(r — 2M) 
-30M^V - 260M^V - 15MV - 15r®) 
and the lowest-order coefficients are: 


(1840M® - 48MV 
(133) 
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0:53 = — 
0:54 = — 
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3a3 

70^3 
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(134) 


= - (^1 - 

RS2 _ ^ 
yrr ^ 5 


2Mr- f3 


= s, 


.,RS2 


2 I 2 I P 2 ■ 2 


= [r^ + a" + ^ 

j^g 2 a{2Mr — (3) sin^ 6 

9 t((> 5 


V sin^ 9 I sin^ 9 


(139) 


where 


A = A- 


(140) 


and where /3 is the tidal charge which can be positive or 
negative. 

This metric is identical to the Kerr-Newman metric 
where the tidal charge is simply the square of the elec¬ 
tric charge Qei- Therefore, this metric could be trivially 
included in the class of metrics designed in this paper 
by starting with the Kerr-Newman metric instead of the 
Kerr metric in Eq. and replacing Qgj —>■ /3. 


VII. DISCUSSION 


C. Chern-Simons Gravity 


Slowly rotating black holes in dynamical Chern-Simons 
gravity were analyzed in Ref. In these solutions, 


In this paper, I designed a Kerr-like black hole metric 
which is regular, admits three independent, exact con¬ 
stants of motion, and depends on four deviations func¬ 
tions in a nonlinear manner. This metric contains the 
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Kerr metric as the special case when all deviations van¬ 
ish. This metric does not derive from any particular 
set of field equations, but, instead, can be used as a 
framework for model-independent strong-held tests of the 
no-hair theorem with observations of black holes in the 
electromagnetic spectrum. I showed that the event hori¬ 
zon of this metric is identical to the event horizon of 
the Kerr metric in Boyer-Lindquist coordinates and that 
the Killing horizon coincides with the event horizon. I 
also determined the range of the deviation parameters 
for which the event horizon and the exterior domain are 
free of pathologies such as singularities or closed timelike 
curves. 

When expanded to linear order in the deviation param¬ 
eters, this metric can be related to the Kerr-like metric 
of Ref. M- I found an explicit general mapping between 
the new metric and the metric of Ref. |14l | in the form 
found in Ref. [l^ , which holds at least up to the fifth or¬ 
der in an expansion of the metric elements in 1/r. I like¬ 
wise mapped my new metric to known four-dimensional, 
analytic black hole solutions in modified theories of grav¬ 
ity. 

I showed that the equations of motion for a test par¬ 
ticle on a geodesic orbit in this metric can be written in 
first-order form and found the Carter-like constant which 
is the third constant of motion along the orbit of the par¬ 
ticle in addition to its energy and axial angular momen¬ 
tum. For particles on circular equatorial orbits, I derived 
expressions for its energy, axial angular momentum, and 
dynamical frequencies and I calculated the location of the 
ISCO. I showed that these quantities depend significantly 
on the deviation parameters. 

These properties make this metric a well-suited frame¬ 
work for strong-held tests of the no-hair theorem in the 
electromagnetic spectrum. Thanks to the existence of an 
event horizon and the regularity of the exterior domain, 
both geometrically thin and thick accretion hows can be 
properly modeled (see the discussion in Refs. [T^. Is^l. 
The hrst-order form of the equations of motion increases 
the speed and precision of ray-tracing codes that are used 
to model and predict observational signatures of non- 
Kerr black holes, because these otherwise have to solve 
the second-order geodesic equations (see Refs. [H, H^). 
In addition, for tests of the no-hair theorem in the 
gravitational-wave spectrum, the mapping of the metric 
linearized in the deviations parameters establishes the 
connection to the metric of Ref. [l3| , for which approxi¬ 
mate EMRI were constructed in Ref. [l^. 

Potential deviations from the Kerr metric should be 
observable with several techniques in the electromagnetic 
spectrum including the continuum-fitting and iron line 
methods. See Refs. for comprehensive reviews on 
these methods. These two methods directly measure the 
location of the ISCO. Since the ISCO does not depend on 
the deviation function (r), this type of deviation would 
be difficult to detect with these techniques. Such partic¬ 
ular deviations, however, may be detected with very-long 


baseline imagin g ob servations of supermassive black holes 
(c.f., e.g.. Ref. [l3|) or of quasi-periodic variability (c.f., 
e.g., Ref. dl]). 

In order to facilitate fully relativistic magnetohydrody¬ 
namic simulations of accretion flows in the new metric, I 
constructed a transformation to Kerr-Schild-like coordi¬ 
nates, which properly removes the coordinate singularity 
at the event horizon in Boyer-Lindquist-like coordinates. 
I demonstrated that the lapse is positive across the event 
horizon and that the amount of frame-dragging near the 
black hole depends strongly on the deviation parameters, 
which can greatly amplify the amount of frame-dragging 
induced by the spin alone. Such accretion flow simu¬ 
lations are carried out in a stationary black hole back¬ 
ground where the dynamical properties of the gravity 
theory do not need to be known. Instead, only the cou¬ 
pling of matter to electromagnetic fields in the modified 
theory needs to be specified. A first step is to assume that 
such interactions are governed by the laws of Maxwell 
electrodynamics. 

Refs. [ 4 ^ 113 implemented a general-relativistic mag¬ 
netohydrodynamic code that is based on a stationary 
black hole metric and several technical assumptions 
which, among others, allow for the treatment of the 
plasma as a perfect fluid. These assumptions also re¬ 
quire that the exterior domain endowed with the metric 
is globally hyperbolic. While this property is not proven 
here, it seems plausible that global hyperbolicity holds 
in the Kerr-like metric designed in this paper, because 
the exterior domain is regular, and, thus, it should be 
possible to trace any flow particle outside of the event 
horizon uniquely to future infinity (or the horizon) for 
any set of initial conditions. Simulations of this nature 
should, therefore, also be feasible in this metric and will 
be explored in detail in a future paper. 
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Appendix A: General Form of the MGBK Metric 

Here I write explicitly the general form of the metric 
of Ref. , where a slightly different notation was used. 
In the notation I use in this paper, the metric is given by 
the elements 


21 


^MGBK _ _____£A)MGBK 


a 2 a^r{r^ + a‘^)Asm9 cos9 


2M Pi 


2a^r'^Asiv?9 jp{r'^ 

P, -+- 

I]A2sin2 0d7i 


dr 

a^)A 


SPi 


SPi 

a A sin^ 9 P 3 
“3 ——FF—Tr73 


L MGBK 
^r6> 


(r^ + a^)/3^A, 


M 


Pi 


2 AP 4 


a A^ {Yi+ 20 ^Mr sir? 9) sir? 9 d'y^ a? A^(E — 4Mr) sin^ 0 ^74 


2M Pi 

tMGBK _ 1 - 7 - 1 

^rr ^3 ? 


dr 2M 


‘■y>y> 


(r^ + or) 


212 




A 2(r2 + 
a 2 a 


a^) 


Pi 


^MGBK 

^t(p 


, 2a2sin6»cos6»^„<,BK , 

— ~^ee B- 


dr 
dr^ 


'-re 


+ 2 - 


90 


dr 2M 


A ^ 2A2 sin^ 6 » 

H- 2 ^3-73 

a 

2 r ^ ^ 2 r„ 

—P - 2r7i + —03 , 


Pi 


2 A^ 
^-74 : 


dr 


8 aM sin 9 cos 9 P 5 


S4 


7e 


4aMr(r^ + a^) sin 6 * cos 6 * 2a^ sin^ 9 Pq 


2^ P 7 7_MGBK 


S Pi 

P, 




2 sin 9 Pio 
H-7771- 77-73 


Pi 

4aMrsin^0Pi5 
^ S2 7^ 
16aM sin^ 9 Pn 


E3 

4aMrsin^0P8 
P-7771-—71 


dr E2 

4aM r Pg ^ 
"^03 


Pi 


lmgbk 


-74 ■ 


E2 YP. P^ 

2a sin^ 9 P 14 ^74 aA sin^ 9 d^ 7 i A sin^ 9 


E2 P^ £2 p^ 

2a P 12 ^71 2 sin^ 9 P 13 ^73 


E2 Pi dr E 

aA(E — AMr) sin^ 9 d?pi 


dr"^ 


E2 


E2 

(E + 2a2Mrsin2 6») 


Pi dr 
(P 73 
dr^ 


E2 


dr2 


(Al) 


where 


03 = 03 ( 0 ), 

7i = 7i(r), i = l,3,4, 


are arbitrary functions of the polar angle 9 and the radius 
r, respectively, and 


(A2) 

(A3) 


I = dr 


2a? sin 9 cos 9 , 


.MG 


2 r 71 + E — 
dr 


(A4) 

The function p is given in Eq. (jll9l) and the functions Pj , 
j = 1 — 15, are polynomials in r and cos 0, given explicitly 
in Appendix A of Ref. H- 
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